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SUMMARY 


An interactive graphical proximity operations planning system has been devel- 
oped, which allows on-site design of efficient, complex, multiburn maneuvers in a 
dynamic multispacecraft environment. Maneuvering takes place in and out of the 
orbital plane. The difficulty in planning such missions results from the unusual 
and counterintuitive character of orbital dynamics and complex time-varying opera- 
tional constraints. This difficulty is greatly overcome by visualizing the relative 
trajectories and the relevant constraints in an easily interpretable graphical 
format, which provides the operator with immediate feedback on design actions. The 
display shows a perspective bird's-eye view of a Space Station and co-orbiting 
spacecraft on the background of the Station's orbital plane. The operator has 
control over two modes of operation: (1) a viewing system mode, which enables the 

exploration of the spatial situation about the Space Station and thus the ability to 
choose and zoom in on areas of interest; and (2) a trajectory design mode, which 
allows the interactive "editing" of a series of way points and maneuvering burns to 
obtain a trajectory that complies with all operational constraints. A first version 
of this display has been completed. An experimental program is planned in which 
operators will carry out a series of design missions which vary in complexity and 
constraints. 


INTRODUCTION 


The future Space Station environment will include a variety of spacecraft 
co-orbiting in close vicinity. Usually, these spacecraft will be placed in a stable 
location with respect to the Space Station; i.e., they will share the same circular 
orbit. However, some missions will require repositioning of or transfers to and 
from these spacecraft. In this case complex maneuvers are anticipated which involve 
a variety of spacecraft not necessarily located at stable locations, and thus there 
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will be relative motion between them. (See appendix for an expanded treatment of 
the subject. 

The multivehicle environment poses new requirements which do not exist in 
conventional mission scenarios. The conventional scenarios involve proximity opera- 
tions between only two vehicles. In these two-spacecraft missions, the scenario is 
optimized through considerable advanced planning and executed at the time of the 
actual mission. However, since the set of possible scenarios in a multivehicle 
environment is virtually unlimited, the future Space Station environment will create 
scenarios that might not have been precomputed and that will have to be planned and 
executed on-site. This will require an on-site planning tool which allows, through 
a fast interactive process, the creation of a fuel-efficient maneuver which meets 
all constraints set by safety rules. 

The difficulties encountered in planning and carrying out orbital maneuvers 
originate from several causes. The first is the counterintuitive character of 
orbital motions as experienced in a relative reference frame. The orbital motions 
are expressed in a coordinate frame attached to the Space Station and represent 
relative rather than absolute motions. It would be intuitively assumed that a 
thrust in the forward direction, i.e., in the direction of the orbital velocity 
vector, would result in a straightforward motion. However, after several minutes, 
orbital mechanics forces will dominate the motion pattern and move the spacecraft 
upward, i.e., to a higher orbit. This will result in a backward relative motion, 
since objects in a higher orbit move more slowly. Thus, a forward thrust has an 
effect opposite from that intended. 

A second cause of the difficulty is the different and unusual way in which 
orbital maneuvering control forces are applied. In atmospheric flight, control 
forces are applied continuously to correct for randomly appearing atmospheric dis- 
turbances, or to compensate for atmospheric drag. In contrast, spaceflight in the 
absence of atmospheric disturbances has a near-deterministic character. Therefore, 
spaceflight is mainly unpowered along a section of an orbit with certain character- 
istics. By applying impulse-type maneuvering forces at a given way point, the 
characteristics of the orbit are altered. After application of the maneuvering 
force, the spacecraft will coast along on the revised orbit until the next way point 
is reached. 

Third, multivehicle orbital missions are subject to stringent safety con- 
straints, such as clearance from existing structures, allowable approach velocities, 
angles of departure and arrival, and maneuver ing-burn restrictions due to plume 
impingement. Design of a fuel-efficient trajectory which satisfies these con- 
straints is a nontrivial task. 

It is clear that visualization of the relative trajectories and control forces 
in an easily interpretable graphical format will greatly improve the feel for orbi- 
tal motions and control forces and will provide direct feedback of the operator's 
control actions. Furthermore, visualization of the constraints in a symbolic graph- 
ical format will enable an interactive graphical trajectory design in which, in each 
iteration step, the design is modified until all constraints are satisfied. 
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BACKGROUND ON ORBITAL MOTION 


Definition of Coordinate Systems and Orbit ^ 

The trajectory of a spacecraft orbiting a given planet is determined by a set 
of six orbital parameters. In figures 1(a) and 1(b) the orbit is visualized. The 
orientation of the plane of the orbit in the inertial system x 1 y 1 z 1 is deter- 
mined by the inclination i and the longitude fi of the ascending node AN (see 
fig. 1(a)). In the orbital plane the shape of the orbit can be either circular, 
elliptic, parabolic, or hyperbolic, depending on the total energy. In our case all 
orbits are either circular or elliptic. Figure 1(b) shows the shape of the orbit in 
the orbital plane. The ellipse is located with one of its foci at the Earth 

center of mass. The orientation of the ellipse is determined by the angle of the 
line of apsides w, which is the angle between the ellipse major axis and the line 
o - AN (which is the intersection of the orbital plane with the horizontal x 1 y 1 
plane). The ellipse is characterized by the semimajor axis a and the eccentric- 
ity e. For a circular orbit the eccentricity is zero and a will be the radius of 
the circle. For the Earth orbit the point closest to Earth is the perigee PE and 
the point farthest away is the apogee AP. Finally, the position of the spacecraft 
in orbit is fully determined by the time of perigee passage t = T, or equivalently, 
by the angle of the true anomaly v. 

In our case all orbital maneuvering takes place about the Space Station. The 
station is stabilized in circular orbit with inclination i = 28.5° with respect to 
the equatorial plane, an altitude of h = 480 km above the Earth surface, and an 
orbital velocity of about V Q = 7,623 m/sec. Trajectory design can be greatly 
simplified by expressing the positions and velocities of co-orbiting spacecraft 
relative to a Space Station-based coordinate system. This system x°y°z° has its 
origin at the center of mass of the station, is oriented with the x°y° plane 
locally level with the surface of the Earth, with the x°-axis in the direction of 
the station's orbital velocity vector, and with the z°-axis pointing toward the 
center of the Earth. Thus, the x°z° plane constitutes the orbital plane. 

The orbital parameter set of a co-orbiting spacecraft can be fully determined 
by measuring the momentary relative position and velocity with respect to the Space 
Station system. These measurements can be performed with great accuracy with Space 
Station-based sensors. The equations for computing the orbital parameters of a 
co-orbiting spacecraft from its measured position and motion relative to the Space 
Station are derived from references 1-3 and are given in the appendix. 


^The coordinate system used here, as well as that in the source code which it 
describes, is the same used in references 1 and 2. It differs from that used in 
many NASA documents in that +R-bar is away from the Earth rather than toward it. 
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Definition of V-Bar, R-Bar, Stationary and Nonstationary Locations 


The V-bar is defined as a section of the circular orbit, followed by the center 
of mass of the Space Station, and the R-bar is defined as the radial line moving 
outward from the Earth center through either the spacecraft or the Space Station 
(see fig. 2). The most appropriate way to determine the in-orbital-plane position 
of a spacecraft relative to the Space Station would be in terms of s and r, where 
s is the distance measured along the circular V-bar between the Space Station and 
the spacecraft's R-bar, and r, the distance of the spacecraft below the V-bar, 
measured along the spacecraft's R-bar (see fig. 2). The velocity vector of the 
spacecraft relative to the Space Station is denoted by v, and its components in the 
V-bar and R-bar direction are denoted by s and r, respectively. Equivalently, the 
vector v can be expressed as the magnitude v and its angle with the V-bar o 
(see fig. 2). Since the range of operation about the Space Station is less than 
several kilometers, the V-bar can be considered to be straight and equivalent to 
the x° axis. Thus x° » s and z° « -r. The out-of-plane displacement is 
y° = w. The relative position is denoted concisely by the vector x. 

A spacecraft, located on the V-bar and with zero velocity relative to the Space 
Station, will, on the assumption of identical drag characteristics, be stationary 
with respect to the Space Station, since it has the same orbital radius R Q and 
period as the Space Station. A spacecraft with radial separation, i.e., above or 
below the V-bar, will not be stationary relative to the Space Station. A spacecraft 
in circular orbit above the V-bar, with radius R = R Q + r, will tend to fall 
behind, i.e., move backward with respect to the orbital flight direction. In con- 
trast, a spacecraft in circular orbit below the V-bar will tend to pull ahead. This 
relative motion is caused by a difference in mean angular orbital rate, or mean 
motion, between the spacecraft and the Space Station. It is well known from 
Kepler's third law that the higher the orbit, the slower the mean motion. If the 
mean motion of the spacecraft in rad/sec is denoted by n, and of the Space Station 
by n Q , the difference in mean motion An is approximated by 



It is seen from equation (1) that after one Space Station orbit, the spacecraft has 
fallen behind or pulled ahead by the distance 3nr meters, where r is the distance 
above or below the V-bar in meters. 

Even when the relative velocity of a spacecraft with radial separation is 
initially zero, it will not be stationary relative to the Space Station. The shapes 
of the orbits for a location above and below the V-bar with an initial radial sepa- 
ration r and an initial zero relative velocity are shown in figure 3(a), and the 
corresponding relative trajectories are shown in figure 3(b). The spacecraft above 
the V-bar is initially at perigee and the spacecraft below the V-bar is at apogee. 
The spacecraft above the V-bar moves initially radially outward away from the Earth, 
and after that starts to fall behind. In contrast, the spacecraft below the V-bar 
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will initially move radially inward toward Earth and after that pull ahead. The 
eccentricity of the orbit and the difference in mean motion between spacecraft and 
Space Station are given by 



( 2 ) 


An » 


' 6 fc) n ° 


and the relative motion starting from t = t Q is approximated by 

s(t) = s(t Q ) - 6r[n Q (t - t Q ) + sin n(t - t Q )] 

r(t) = r(t ) + 3r[1 - cos n(t - t )] 
o o 


(3) 


(4a) 

(4b) 


Equations (4a) and (4b) are obtained by solving the general equations for relative 
motion given in the appendix for the particular case in which r(t Q ) * 0 and s(t Q ) 
and r(t Q ) are zero. It is clear from equations (4a) and (4b) that after one Space 
Station orbit the spacecraft is again back on its original radius R, but it has 
fallen behind or pulled ahead by the distance 12nr meters where r is the distance 
above or below the V-bar in meters. The maximum deviation in r-direction from 
r(t Q ) occurs at one-half orbit and is 6r meters. 


Typical Maneuvers — V-Bar and R-Bar Burn 

Consider a spacecraft located at the V-bar, and thus at a stationary position 
relative to the Space Station. A small maneuvering burn in the direction of the 
positive R-bar will cause a small component v which will result in a small change 
in the direction of the orbital velocity vector. This will alter the parameters of 
the orbit. The orbit will become elliptical and the spacecraft will be at apogee 
after one quarter of an orbital period. In figures 4(a) and 4(b) the shape of the 
orbit and the corresponding relative motion trajectory is shown. The relative 
trajectory has a closed elliptical shape and after one orbit the spacecraft will 
return to its original location. The reason for this is that the radial burn did 
not significantly alter the magnitude of 
motion did not change. The eccentricity, 
motion for an R-bar burn are given by 


e = 


An « 


V Q and thus the total energy and mean 
difference in mean motions and relative 


v_ 

V 


(5) 


( 6 ) 
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(7a) 


s(t) * s(t ) - 2R 
o o 



[1 - cos n (t - t )] 
o o 


r(t) * R 


ft) 


sm 


n 0 ( fc - 


V 


(7b) 


Equations (7a) and (7b) are obtained by solving the general equations in the appen- 
dix for the case r(t Q ) = 0, s(t Q ) = 0, and r(t Q ) = v * 0. Since n Q = V Q /R 0 , it 
is seen from equations (7a) and (7b) that for a positive R-bar burn of 1 m/sec, 
after half the orbital period of the Space Station, the spacecraft is again at the 
V-bar, but has fallen behind the distance 4/n » 3,598 m. The maximum deviation 
from the V-bar occurs at one-fourth and three-fourths of an orbital period of the 
Space Station and is 1/n a 899 m. 

In contrast to the R-bar burn, a maneuvering burn along the V-bar will alter 
the magnitude of V Q by the amount v and will therefore alter the total energy. 
Figures 5(a) and 5(b) show the shape of the orbit and the corresponding relative 
motion trajectory. For a burn in the positive V-bar direction, the spacecraft will 
initially move forward, but later on it will gain altitude and fall behind. The 
opposite is the case for a burn in negative V-bar direction. Here the spacecraft 
will initially move backward, but later on drop altitude and pull ahead. For a 
positive burn the spacecraft is initially at perigee and for a negative burn at 
apogee. The eccentricity, difference in mean motions, and relative motion for a 
V-bar burn are given by 


e = 




( 8 ) 


An » 



S ( t ) S s(t ) + R 
o o 


ft) '•». 


(t - t ) + 4 sin n(t - t ) ] 
o o 


r(t) * 2R 


ft) 


[1 - cos(t - t )] 
o 


(9) 

(10a) 

(10b) 


Equations (10a) and (10b) are obtained by solving the general equations in the 

appendix for the case r(t Q ) = 0, s(t Q ) = v * 0, and r(t Q ) = 0. Equations (10a) 

and (10b) show that for a positive V-bar burn of 1 m/sec, after one orbital period 

of the Space Station, the spacecraft has fallen behind the distance 

6ir/n » 16,957 m and the maximum distance from the V-bar is 4/n a 3,598 m. 
o o 
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Limitations of Conventional Techniques 


It is clear from the above examples that orbital motion can be complex and 
highly counterintuitive. A burn toward the target might have the opposite result 
than intended. Motion is particularly complex and difficult to visualize for a 
combined burn at a nonstationary location. It is therefore very useful to graphi- 
cally visualize the relative motion trajectories. This is true in particular for 
on-site interactive trajectory design in which immediate graphical feedback on the 
results of maneuvering burns is essential. 

The presently used conventional techniques are well established and rely in 
most cases on visual contact and the use of a V-bar or R-bar reference in a Crewman 
Optical Alignment Sight (COAS) (ref. 4). In a V-bar approach toward a target in 
positive V-bar direction, the initial burn is made in a direction slightly depressed 
downward with respect to the V-bar. After a short while, the spacecraft will ascend 
and cross the V-bar. At the V-bar crossing, a small downward R-burn is initiated 
which again depresses the spacecraft below the V-bar. This process is repeated 
several times. The spacecraft thus proceeds along the V-bar in small hops until the 
target is reached. A similar approach is possible along the R-bar. However, these 
scenarios are highly restricted in their complexity and are not likely to satisfy 
the complex constraints of the multivehicle environment. 


DESCRIPTION OF THE TECHNIQUE 


Purpose of Orbital Planning System 

The purpose of the interactive orbital planning system is to enable the opera- 
tor to design an efficient, complex, multiburn maneuver, subject to the stringent 
safety constraints of the future dense traffic environment of the Space Station, 
which enables a chaser to rendevous with a target spacecraft in a given timespan. 

The constraints include clearances from structures, relative velocities between 
spacecraft, angles of departure and arrival, approach velocity, and plume impinge- 
ment. Because of the complexity and counterintuitiveness of orbital motion, and the 
demands to satisfy strict safety rules and constraints, fuel-efficient trajectory 
design will be a complex and difficult task. The basic idea underlying the system 
is to present the maneuver as well as the relevant constraints in an easily inter- 
pretable graphical format. This format provides the operator with immediate feed- 
back on the results of design actions, and allows for close interaction with the 
system. In an iterative process, the operator will keep changing the design until 
all constraints are met. The methods for enabling interactive trajectory design and 
visualization of constraints are discussed in detail hereafter. 
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Illustration of a Three-Burn Maneuver 


An illustrative example of a three-burn maneuver is shown schematically in 
figure 6. The trajectory originates from relative position A at time t = t Q , and 
is composed of two way points B and C which specify the location in Space Station 
coordinates through which the chaser spacecraft will pass at a given time. At a way 
point the orbital maneuvering system or other reaction control system is activated, 
creating a thrust vector of given magnitude for a given duration, in a given direc- 
tion in orbital plane or out of the orbital plane. The duration of the burn is 
considered very short in comparison with the total duration of the mission. In the 
orbital dynamics computations this means that a maneuvering burn can be considered 
as a velocity impulse which alters the direction and magnitude of the instantaneous 
orbital velocity vector of the spacecraft. 

Since the initial location A is not necessarily a stationary point, the 
magnitude and direction of the relative velocity of the chaser at point A is 
determined by the parameters of its orbit. If no maneuvering burn is to be ini- 
tiated at t = t Q , the chaser will continue to follow the relative trajectory 1, 
subject to the parameters of its original orbit (see dotted line in fig. 6). How- 
ever, a maneuvering burn at t = t Q will alter the original orbit; thus the chaser 
will follow the relative trajectory 2, subject to the parameters of the new orbit. 

In figure 6, v-j and v 2 indicate the relative velocity vector of the chaser 
just before and after the maneuvering burn, where v.| and v 2 are tangential to the 
relative trajectories 1 and 2, respectively. The vector difference between v^ and 
v 2 , ^ is the velocity change initiated by the burn, and corresponds to the direc- 
tion and magnitude or duration at which the orbital maneuvering system is acti- 
vated. Likewise, at way point B the burn v b alters the orbit to orbit 3. 

Location C is the terminal way point and in this case is the location where 
the target will arrive at t = t^. Since the target has an orbit of its own, 
orbit 4, it will have a terminal velocity at t = t^. The relative velocity between 
target and chaser is the vector difference between v^ and v^, v c . This vector 
determines the retro-burn that is needed at the target location to bring the rela- 
tive velocity between chaser and target to the minimum required for the docking 
operation. 


Inverse Method of Solving Orbital Motion 

Interactive trajectory design demands that the operator is given free control 
over the positioning of way points. However, the input variables of the equations 
of orbital motion given in the Appendix are the magnitude and direction of the burn, 
rather than the position of way points. Therefore an inverse method is required to 
compute the values of a burn necessary to arrive at a given way point positioned by 
the operator. This method is outlined here. 

The equations in the appendix show how the orbital parameters of a co-orbiting 
spacecraft can be computed from its momentary position and velocities, relative to 
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the Space Station. Thus, for a given initial relative position A with x(t 0 ) anc * 
an initial relative velocity v(t Q ) at time t = t Q , the relative position and 
velocities of a way point at time t = t^ can be computed. However, a maneuvering 
burn at t = t Q will cause a change in the direction and magnitude of the relative 
velocity vector v(t Q ). As a result, the position of the way point at time t = t^, 
x(t^) will change as well. 

Consider v a and a a to be the magnitude and direction of the velocity change 
due to the maneuvering burn. Then the relative position and velocity at t = t -j , 

x(t^), will be a complex nonlinear function of v a and a a . Consider now that the 

operator is given direct control over v a and a a by slaving these variables 
directly to the x and y motions of the mouse, respectively. A displacement of the 
mouse in either x or y direction will result in a complex nonlinear motion pattern 
of x(t^). Furthermore, this motion pattern will change with the initial condi- 
tions. This arrangement is highly undesirable in an interactive trajectory design 

process in which the operator must have direct and unconstrained control over the 

positioning of way points. 

It is therefore essential to give the operator direct control over the position 
of way points rather than over the magnitude and direction of the burn. The inverse 
method by which this is accomplished computes the magnitude and direction of the 
burn required to bring the spacecraft from initial location x(t Q ) at t = t Q to 
the way point x(t^) at t = t-j. A Newton-Raphson method has been employed to solve 
this inverse problem. The operator commands the position of a way point by means of 
the x-y motions of the mouse. The algorithm starts with an initial guess of v a 
and a a . These values yield a computed way point which is usually different from the 
commanded one. At each program update the values of v a and a a are adjusted to 
bring the computed way point closer to the commanded one. On the average about 
three or four iterations are required to bring the difference between the computed 
and commanded way point effectively to zero. As the operator moves the commanded 
way point around in the orbital plane, the algorithm tracks the commanded way point 
by continuously making appropriate adjustments in v a and a a . Because of this 
continuous adjustment, the deviation between commanded and computed way point will 
remain relatively small and the Newton-Raphson scheme will operate in the vicinity 
of the optimum. The advantage of the Newton-Raphson scheme is that convergence with 
this second-order technique is the best in the near vicinity of the optimum. Since 
the program update rate is about 15 Hz, convergence is very fast and the computed 
way point is virtually indistinguishable from the commanded one. 


Concept of the Active Way Point 

Although a trajectory may be composed of several way points, only one way point 
at a time, the active way point, is controlled by the operator. While the position 
and time-of -arrival of the active way point can be varied, the position and time-of- 
arrival of all other way points remain unchanged. However, variations in the active 
way point will cause changes in the trajectory sections and way-point maneuvering 
burns just preceding and just following the active way point. The on-line solution 
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of the inverse algorithm enables these changes to be visualized almost instanta- 
neously and provides the operator with on-line feedback on design actions. 

Although impingement constraints and approach velocity limits exist for all way 
points, it is useful to limit the computation and display of these constraints to 
the active way point only. This arrangement simplifies and speeds up system update 
computations and minimizes the symbology shown on the display. The justification 
for this is that the operator's attention is mainly allocated to the active way 
point and its near vicinity. In a subsequent design iteration the operator may 
shift the activation to a different way point and again verify whether all con- 
straints are met. 

Since impingement constraints and approach velocity limits mainly relate to the 
target craft, it is useful to visualize the position of the target on the target 
trajectory, corresponding to the time of arrival at the active way point. Both the 
active way point, and the target position at active way-point time should be clearly 
distinguishable from other way points by conspicuous markings, highlighting, or 
blinking. 


Process of Way-Point Editing 

The trajectory-design process involves changes in existing way points, addition 
of new points, or deletion of existing undesired points. An illustrative example of 
this way-point editing process is shown in figure 7. In the program the way points 
are managed by a way-point stack, which includes an up-to-date sequential list of 
the position x, the time-of-arrival t, and the relative velocity v just after 
initiating the burn, of all way points. 

Figure 7(a) shows two way points, the initial point Xg and the terminal 
point x^ . The initial way point is defined by the initial conditions of the situa- 
tion and cannot be activated or changed by the operator. The terminal way point 
x^ is thus the active way point which can be changed. The corresponding way-point 
stack is shown on the right. The active way-point box is drawn in bold. The rela- 
tive velocity stack shows only the velocity Vg, which is the relative required 
velocity just after the burn at way point 0, computed by the inverse algorithm, and 

required to reach point Xi at time fc i* 

Figure 7(b) shows the addition of a new way point. The new way point is added 

half-way on the trajectory section just preceding the active way point. Thus its 

time-of-arrival is chosen to be t = 0.5(t i + t^), where i is in this case 1 and 
relates to the stack before modification. The new position x^ and relative veloc- 
ity v 1 are computed by the forward method given in the appendix, by computing the 
orbital position at the new time t using the existing orbital parameters pre- 
viously computed with Xg, Vg, and t Q . The newly computed way-point position, time, 
and relative velocity are inserted between points 0 and 1 of the stack before modi- 
fication, and the new way point is chosen to be the active one. The dotted lines in 
figure 7 indicate variables which are transferred without modification and the 
encircled variables are the newly computed ones. It is important to note that since 
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the relative velocities Vg and v^ are matched to the required way points and 
x 2 » respectively, the inverse algorithm does not need to make any adjustments. 

Figure 7(c) shows the results of changes in the newly created way point on the 
way-point stack. Since x 1 and t 1 are varied, the relative velocity at way-point 0, 
Vq, will be readjusted by the inverse algorithm and likewise the relative velocity 

-1 * 

Figure 7(d) shows the creation of an additional new way point. Since the 

active way point prior to the addition was point 1 , the new point is added half-way 

between points 0 and 1 and its position and relative velocity are computed with the 

forward method. The new values are inserted between points 0 and 1 of the stack 

before modification, and the new way point is again set to be the active one. 

In figure 7(e) way point 2 is activated. Apart from the shift in active way 
point, the stack remains unchanged. The dotted line shows the direct-path section 
between point 1 and point 3 without the intermediate burn at point 2. Deletion of 
way point 2 will remove this point from the stack and, after that, close the gap 
(see fig. 7(f)). However, v.j has to be readjusted to fit the new direct-path 
section. Starting from the old incorrect value of v^ , one can make the adjustment 
iteratively and on-line by the inverse algorithm. 


Constraints During Operations 

The multispacecraft environment will require strict safety rules regarding 
clearance from existing structures. Thus, spatial envelopes can be defined through 
which the spacecraft is not allowed to pass. These spatial constraints can be 
visualized on the display. The operator must be able to make a clear judgment as to 
whether the planned trajectory clears the spatial constraint. In case the trajec- 
tory does not clear the constraint, the operator must be able to decide whether to 
avoid the constraint through an in-plane maneuver or an out-of-plane maneuver. 
However, the operator is not always able to make these judgments on the basis of 
one-perspective aerial view or one-perspective projection. In this research, a 
graphical enhancement is used in which the spatial constraint is unambiguously 
presented in a time-axis display format. This format and its advantages are dis- 
cussed in the "Description of the Display" section of this paper. 

Restrictions on angles of departure and arrival may originate from structural 
constraints at the departure gate or from the orientation of the docking gate or 
grapple device at the target craft. Limits for the allowable angles of departure or 
arrival can be visualized on the display. In addition, the terminal approach veloc- 
ity at the target might be limited by the characteristics of the grapple mechanism 
or the docking procedure. Limits for the allowable terminal approach velocity can 
be visualized as well. 

Way point maneuvering burns are subject to plume impingement constraints. Hot 
exhaust gases of the orbital maneuvering systems may damage the reflecting surfaces 
of sensitive optical equipment such as telescopes or infrared sensors. Maneuvering 


burns toward these pieces of equipment are restricted in direction and magnitude; 
limits for the allowable direction and magnitude are a function of the distance to 
the equipment and plume characteristics. These limits can be visualized on the 
display. 

Flight safety requires that the relative velocity between spacecraft is subject 
to approach velocity limits. In conventional docking procedures these limits were 
proportional to the range (refs. 5 and 6). A commonly used rule of thumb is to 
limit the relative approach velocity to 0.1 percent of the range. This conventional 
rule is quite conservative and originates from visual procedures in which large 
safety margins are taken into account to correct for human or system errors. 

Although the future traffic environment will be more complex, and will therefore 
demand larger safety margins, more advanced and reliable measurement and control 
systems will somewhat relax these demands. The effect of these developments on the 
allowable approach velocity limits is at present difficult to predict, and so is the 
margin for human error to be taken into account. 

In this study, the relative approach velocity is defined as the component of 
the relative approach velocity vector between the two spacecraft along their mutual 
line of sight. The limit on this relative approach velocity is a function of the 
range between the spacecraft. This function will depend on the environment, the 
task, and the reliability of measurement and control equipment, and cannot be deter- 
mined at this stage. In this study a simple proportional relation has been 
chosen. The approach velocity limit is visualized on the display as a circle indi- 
cating the minimum range between the two spacecraft allowed for the present approach 
velocity. If the target craft appears within this circle, the approach velocity 
limit has been violated. 


DESCRIPTION OF GRAPHICS SYSTEM AND LAYOUT OF THE DISPLAY AREA 


The system has been implemented on a Silicon Graphics IRIS 2400 Turbo Graphics 
Workstation with 24 bitplanes of display memory and with a 19-in. full-color display 
monitor with a display resolution of 1024 by 768 pixels. The program is named NAVIE 
(the Hebrew word for prophet) after the prophet Elijah, who was characterized as 
providing trustworthy future information. Operator interaction with the system is 
through a two-axis, three-button mouse. 

The layout of the display area is shown in figure 8. The display area has been 
divided into four viewports, main area 1 (750 by 750 pixels) and areas 2, 3, and 4 
(230 by 230 pixels each). Viewports 1, 3, and 4 provide information about the 
spatial situation around the Space Station, trajectories, constraints, and orbital 
maneuvering fuel use; and viewport 2 includes an eight-button function control 
panel . 
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DESCRIPTION OF DISPLAY PROGRAM CONTROL MODES 


The program operates in two basic modes. The first one, the system control 
mode, is engaged when no buttons are pressed. In this mode, motions of the mouse 
are controlling actions in windows 1 and 3 (such as viewpoint motion or way-point 
position). The second mode is devoted to function control and is activated by 
pressing the middle button of the mouse. When this button is held down, a cursor 
appears which can be positioned on one of the eight function buttons in area 2. 
Selected buttons are highlighted and the function described by the button is acti- 
vated only after the middle button is released. 

While being switched from system control mode to function control mode, the 
system control actions are frozen and reactivated after return from the function 
control mode. Note that although the mouse values may have changed in the function 
control interval, the system control status remains unaltered. Likewise, the latest 
cursor position in the function control mode remains unchanged in the system control 
interval, although the mouse may have been moved. This feature is very useful when 
' the same button is activated more than once. For example, this is done when tog- 

gling through a sequence. The toggling in this case can be done simply by pressing 
I and releasing the middle mouse button several times, without the need for reposi- 

tioning the cursor. 

The system control mode has two submodes. The first one, the viewing system 
' mode, relates to the main display which shows a perspective view of the Space Sta- 

tion and its surroundings on the background of the Station's orbital plane. In this 
viewing system mode, the operator is able to explore the spatial situation about the 
Space Station and thus choose a viewpoint location and viewing direction which focus 
and zoom in on the momentary area of interest. The second submode is the trajectory 
design mode, in which way points are selected, moved, added, and deleted in order to 
obtain a multiburn trajectory which complies with the given set of constraints. 


Viewing System Mode 

The geometry of the viewing situation is shown in figure 9. The Space Station- 
based coordinate system is x°y°z° with the x° axis coinciding with the orbital 
velocity vector, and the orbital plane is x°z°. Figure 9 shows the orientation of 
the viewing system x e y e z e relative to the Space Station system. The viewing 
system has its origin at point A, the x e axis coincides with the viewing direc- 
tion, and the image plane is perpendicular to the x e axis with the screen axes y s 
and z s parallel to y e and z e . Point B indicates the intersection of the viewing 
axis with the orbital plane. Although the viewing system position, point A, and 
angular orientation are defined by three displacements and three angles (which can 
all be controlled independently), it is useful to constrain the motion to the fol- 
lowing three types. 
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Tethered Motion - In the first type of motion the viewing system tethers about 
point B, which is kept fixed on the orbital grid, while the distance d between 
points A and B, which is the viewing range to point B, is kept constant. The 
tethered motion is controlled by the angles 4> and 0. The viewing axis x e and the 
axis y e are located at all times in the plane P which passes through the point 
B and rotates about the line CC', which is parallel to the x° axis (the 
V-bar). The line BE is also located in the plane P, and is perpendicular to the 
line CC'. i|» is the angle between the y° axis and the line BE, and 0 is the 
angle between BE and the x e axis. Thus, the angles i|> and 0 control the obliq- 
uity of viewing along the orbital plane in z° and x° directions, respectively. 

This tethered type of motion is very useful for the following reasons: (1) While 
the area of interest remains in the center of the display, it allows one to explore 
other possible areas of interest by bringing them into the viewport by changing the 
angles 4* and 0. (2) The line CC' will appear on the screen at all times as a 

horizontal line through the center of the display, and will represent a line paral- 
lel to the V-bar. Thus, while the viewing direction may change, the direction of 
the V-bar is at all times recognizable as the horizontal line passing through the 
center of the display. 

Translational Motion - The second type of motion relates to the position of 
point B in the orbital plane. Here the x°z° coordinates of point B are 
varied, while ip, 0, and d are kept constant. This type of translational motion 
enables the operator to move areas of interest to the center of the display. 

Ranging Motion - In the third type of motion, all parameters are kept constant 
except for the range d. This ranging type of motion is useful after areas of 
interest are located and brought into the center of the display. Ranging-in on the 
area of interest allows this area to be studied in more detail. 

In the viewing system mode, the operator can select either tethered motion, 
translational motion, or ranging motion by pressing the middle button. This one- 
button control is useful since viewing system operations are naturally performed in 
a sequence of three steps. In the first step, areas of interest are selected. In 
the second step the selected area of interest is moved to the center of the display, 
and in the third step the area of interest is ranged-in on to obtain the required 
level of detail. Note that a ranging operation rather than a zoom operation has 
been chosen. This has been done to preserve the operator's awareness of the spatial 
location of the viewpoint. The ranging operation is consistent with the transla- 
tional type of motion of point B, both of which are experienced by the operator as 
flying in three-dimensional space around the Space Station. In a zoom operation the 
viewpoint position remains unchanged and instead, the field of view is changed. 
However, since zooming is experienced as moving toward or away from the viewed 
object, this operation is bound to interfere with the operator's spatial 
orientation. 

In addition to the control over the three types of motion, the viewing system 
mode has one additional control, a viewing system reset. This reset will return the 
viewing situation to a default setting which is defined in a program parameter file 
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prior to the running of the program. The default viewing situation can be chosen to 
view the total area of interest from a more remote location. 


Trajectory Design Mode 

In the trajectory design mode, the operator has control over the seleotion, 
positioning, time of arrival, and addition and deletion of the way points which 
determine the trajectory. Two submodes exist in the trajectory design mode, the 
in-plane design mode, and the out-of-plane design mode. In the in-plane mode the 
mouse controls the x°z° position of way points, while the out-of-plane position 
y° remains unchanged. In the out-of-plane mode the in-plane position remains 
unchanged and the out-of-plane position can be varied. 

The design process starts with an initial configuration of way points. Usually 
there initially are two way points as in the example of the Way-Point Editing sec- 
tion. The terminal point Xi is the active way point. Time of arrival at the 
active way point is set to an initial value within the allowable timespan of the 
mission. The operator can increase or decrease the time of arrival at any active 
way point in preset steps by repeatedly pressing the right or left button, respec- 
tively, of the mouse. In this example a time-step of 10 sec was found to be ade- 
quate. The time of arrival at the terminal way point is limited to the timespan of 
the mission, and the arrival time at intermediate way points is restricted by the 
timespan set by the neighboring points. 

As outlined in the Way-Point Editing section, a convention is chosen in which a 
new way point is added to the trajectory section between the way point preceding the 
active way point, and the active way point itself, half-way on the time scale 
between the two way points. The newly added way point becomes the active one and 
can be moved to any desired location, and its time of arrival can be set to any 
value within the timespan determined by the neighboring way points. However, in 
some cases, it is useful to slide the new way point along the trajectory section 
directly connecting its neighboring way points. The position on this trajectory 
section is then determined by its time of arrival only. In this mode, the "locked- 
on-trajectory" mode, the time of arrival is slaved to the y-motions of the mouse. 

The locked-on-trajectory mode is particularly useful for checking whether 
operational constraints between the spacecraft and the target or other nonstationary 
spacecraft are being violated. As the operator slides the way point along the 
trajectory, the corresponding target position slides along the target trace as well; 
conflicting situations, such as a too close fly-by, will be recognized immediately. 


Geometrical Enhancements — The Time-Axis Format 

The purpose of these enhancements is to resolve ambiguities in the spatial 
situation by processing the spatial information and presenting it in a different 
format. One such format is the time-axis display, which provides unambiguous 
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qualitative and quantitative information about the out-of-plane situation and the 
spatial constraints. 

The basic idea of the time-axis format is demonstrated in figures 10(a)-(c). 
From the perspective view of figure 10(a) alone it cannot be clearly determined 
whether the spatial constraint is violated or how the trajectory should be planned 
to avoid it. The view along the z° axis in figure 10(b) is even less clear, 
because of the curved character of the trajectory. In the time-axis format of 
figure 10(c) the out-of-plane deviation is plotted as a function of the traveled 
time along the path. The spatial constraints are visualized as follows. At each 
point on the traveled-time axis, at the corresponding location on the trajectory, a 
line is placed perpendicular to the orbital plane. Sections of this line which are 
within these constraints are identified and plotted on the time-axis display of 
figure 10(c) as a set of vertical bars. Where the trajectory curve passes through 
these bars, the spatial constraints have been violated. From this display it can 
also be clearly determined whether the constraint should be avoided through an 
in-plane or an out-of-plane maneuver. Depending on the shape of the spatial con- 
straint, an in-plane maneuver might vary the size and location of the constraint 
bars. For the rectangular constraint in figure 10(a) the bars might disappear 
entirely. In other cases, a small out-of-plane maneuver might be sufficient to 
avoid the spatial constraint bars. 

The format of the time-axis display used in the program is shown in fig- 
ure 11. The time axis is marked in quarters of an orbit. The shaded area repre- 
sents the nighttime section of the orbit. This is useful in verifying requirements 
regarding the section of the orbit, i.e., day or night, in which a certain part of 
the mission has to take place. Both the target and the chaser trajectories are 
shown. It should be noted, however, that although the chaser and target share the 
same time axis, they relate to different spatial trajectories. Therefore the spa- 
tial-constraint bars are meaningful for the chaser trajectory only. 

The time-axis display can be brought into the main viewport by a swap between 
viewports 1 and 3. This viewport swap feature is useful when certain parts of the 
time-axis display have to be viewed in more detail and the resolution of viewport 3 
is insufficient. 


Symbolic Enhancements 

Visualization of Departure Constraints - Procedures at the departure gate might 
constrain the relative angle of departure and the magnitude of the departure burn. 
The in-plane constraints at the departure gate are visualized in figure 12. The 
size of the burn vector is made proportional to the burn magnitude. The departure 
constraints are satisfied if the burn vector is within the solid bracketed arc. 

This arc is specified by the arc center angle y 0 » the arc aperture y, and the arc 
radius e. The scale-factor for visualizing the burn velocities in this example is 
500 m length per 1 m/sec burn. Note that maneuvering burns are expressed in terms 
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of a velocity change rather than of a thrust force. The actual duration and thrust 
force of the burn depends on the spacecraft mass and the thruster characteristics. 

In order to keep the display free from unnecessary symbology it is useful to 
present the constraint only when it is close to being violated. If the burn vector 
is within the area enclosed by the dotted line in figure 12, the constraint is not 
drawn. The radius of the dotted arc is 80% of e and the aperture angle is 10° 
smaller than y • 

It should be noted that the situation in figure 12 relates to a stationary 
departure gate. The spacecraft trajectory in this case is aligned with the burn 
vector. For a departure gate which moves with respect to the Space Station system, 
this will not be the case. In this case the burn vector will signify the relative 
direction of departure with respect to the moving gate, rather than with respect to 
the Space Station. But this is just the vector which is subject to the departure 
constraints and not the velocity vector of the trajectory, which is relative to the 
Space Station. Therefore, the symbology is valid for departure from a stationary as 
well as nonstationary gate. 

The out-of-plane constraint at the departure gate is visualized in figure 11. 
The initial out-of-plane component of the burn vector has to be within the impinge- 
ment constraint brackets. If the burn magnitude is less than &0% of the allowed 
maximum value, the constraint is not drawn. 

Visualization of Arrival Constraints - The arrival procedures constrain the 
angle and magnitude of the terminal velocity vector relative to the arrival gate. 

The in-plane constraints at the arrival gate are visualized in figure 13. The size 
of the relative terminal velocity vector is made proportional to the vector 
length. The arrival constraints are satisfied if this vector is within the solid 
arrival arc. This arc is specified by the arc center angle 6 Q , the arc aperture 
6, and the arc radius n. The arrival arc is visualized at all times. The scale- 
factor for visualizing the terminal velocity in this example is 500 m length per 
1 m/sec. 

The out-of-plane limits on the terminal approach velocity are depicted in 
figure 11. The approach velocity has to be within the constraint brackets. If the 
velocity is less than 80^6 of the allowed maximum value, the constraint is not drawn. 

Visualization of Plume Impingement Constraints - Plume impingement constraints 
limit the magnitude and direction of a maneuvering burn. The in-plane impingement 
constraints of a burn given at a way point toward the target are illustrated in 
figure 14. The burn-vector symbol, whose size is proportional to the magnitude of 
the burn, is not allowed to cross the bracketed impingement constraint arc with 
aperture 8 and radius a. The variables 8 and a are a function of the distance 
between way point and target |Ax| = |x T - x|, whose function depends on the charac- 
teristics of plume and target. In this example, 8 is chosen to be constant and 
a proportional to | Ax | . If the burn vector does not cross the dotted bracketed 
arc, the constraint is not drawn. The radius of the dotted arc again is 80% of a 
and the aperture angle 10° larger than 8. The scale factor for visualizing the 
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burn velocities was again chosen in this example to be 500 m length per 1 m/sec 
burn. 


Visualization of Approach Velocity Constraints - The method of visualizing the 
relative approach velocity limit is shown in figure 15. The relative approach 
velocity of the chaser toward the target is given by the vector Av = v - v.p. The 
line-of-sight vector of the chaser toward the target is Ax = “ £• The relative 

approach velocity vector v r is the projection of Av on Ax and is given by 

Ax 

v = 5 (Av 1 Ax) (11) 

~ r i*il 2 ' ' 

where T denotes the transpose, or inner product. The limit on |v r | is a func- 
tion of the distance between chaser and target | Ax | . In this example, a simple 
proportional relationship has been chosen. Thus, for a given approach velocity 
|v r |, the allowable range p can be computed and visualized by a circle centered 
about the chaser's position. The approach velocity constraint is violated when the 
target is located within this circle. The circle is visualized when p is greater 
than 80# of | Ax | . 

Orbital Fuel Use - The orbital fuel use is displayed in viewport 4. The orbital 
fuel is expressed in total velocity change in meters per second rather than by fuel 
use in kilograms. The actually spent fuel mass depends on the spacecraft and the 
thruster characteristics and will be proportional to the total velocity change. A 
fuel dial is shown which indicates the percentage of fuel remaining from the total 
amount allowed for the mission. The remaining fuel is indicated by a yellow sector, 
and fuel use in excess of the allowed amount is indicated by this sector turning 
red. In addition to the fuel dial, the percentage of fuel left and total fuel use 
are displayed numerically. 

Trajectory Time Markers - Time markers are placed at regular intervals along the 
chaser and the target trajectories, indicating the time after the start of the 
maneuver. The time marker is a small bar perpendicular to the trajectory. At the 
end of this bar a number is displayed which indicates the time in minutes after the 
beginning of the maneuver. It should be noted that in the viewing system mode, the 
obliquity and range at which a time marker is located might change. When the num- 
bers are written in the orbital plane, they will change with each viewing system 
change. A worst-case example occurs when the backside of the orbital plane is 
viewed and the numbers appear in mirror script. This can be prevented by writing 
the numbers in the plane of the image rather than in the orbital plane. However, a 
viewing system change will cause the numbers to move with respect to the marker. 
Special care is given to the automatic repositioning of the numbers after a viewing 
system change. The numbers should be placed such that they do not clutter the 
trajectory, and clearly point to the corresponding time marker. 
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Computational Enhancements 


Computation of the relative trajectories is a time-consuming process, which, if 
done at each program update, results in an unacceptably low update rate. This 
results in jerky motions and poor control over the positioning of a way point. This 
problem can be prevented by disabling the trajectory computations and starting them 
only after the operator has stopped moving the active way point. 

In order to avoid disturbing delays, at each program update interval the x and 
y output values of the mouse are compared with the values from the previous step. 

If no change has taken place, a timer is initiated. The trajectory computations are 
only started 0.3 sec after initiating the timer. After the trajectory is computed, 
the computed values are stored and displayed and no further computations will take 
place until the next change in the active way-point position. The 0.3-sec delay is 
essential for ensuring that the operator has completed the positioning process. 

Often small corrections are made after the way point has been moved but these cor- 
rections do not interfere with the program's update rate, because of inhibition of 
the trajectory calculation. Sometimes subsequent changes are made after the opera- 
tor has reviewed the position. These changes are seldom made earlier than 0.5 sec 
after the last change, and this is after the trajectory has been recomputed. 

It should be noted that although the trajectory computations are subject to 
delay, this is not the case with the computation of variables which relate to the 
way points themselves, such as maneuver ing-burn vectors, relative-velocity vectors 
and operational constraints. The computation of these variables is less time- 
consuming and is done at each program-update interval. Continuous update of these 
variables is essential in order to give the operator immediate feedback of the 
effect of a certain design action on maneuvering burns or approach velocities. 

A second enhancement for minimizing the computational load of the trajectory 
computations is a variable time-step size algorithm. The trajectory points are 
computed at fixed intervals of time rather than at fixed distance intervals along 
the trajectory. Since the trajectory length might vary, the distance between the 
trajectory points might vary as well. At each program iteration, the total trajec- 
tory length is computed and divided by the total trajectory time to obtain an aver- 
aged trajectory velocity. The time step is computed by dividing a preset required 
averaged distance between trajectory points by the averaged trajectory velocity. 


CONCLUDING REMARKS 


The proposed interactive orbital planning system should be seen as a prelimi- 
nary step in determining the display format of graphical trajectory design, which 
will be useful in the dense Space Station environment. The examples shown here deal 
with the most general situation, which involves departures from, and arrival at, 
nonstationary locations. However, most of the co-orbiting spacecraft are likely to 
be placed on the V-bar, and thus, at stationary positions. Missions with spacecraft 
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at nonstationary positions and substantial out-of -plane motion thus represent a 
worst-case situation, and as shown here are chosen to demonstrate the capabilities 
of interactive graphical trajectory design, rather than to represent the average 
type of maneuver to be executed at the Station. 

Likewise, it is difficult to predict whether the constraints used here will be 
relevant and realistic in the future Space Station environment. They predict in a 
broad sense the types of restrictions which are expected in the multivehicle envi- 
ronment; e.g., limitations on approach rates, plume impingement, and clearance from 
structures. It is also likely that the future environment will pose different 
constraints that might originate from the specific character of a mission, such as a 
specific scenario in which a telescope or manufacturing platform is approached and 
serviced. 

A third restriction of the display relates to the way the orbital maneuvering 
system is activated. Only pure-impulse maneuvering burns are considered, in which 
the duration of the burn is negligible with respect to the duration of the mission, 
and in which these burns cause major changes in the relative trajectories. Station- 
keeping or fly-by missions, however, require a more sustained type of activation, 
such as periodic small burns with intervals of several seconds over a time span of 
several minutes. A more distributed way of activating the orbital maneuvering 
system can be introduced in which the operator has control over the frequency and 
time span of the activation. Ways can be found to enable this type of control to be 
activated and visualized. 

A last restriction relates to the way the spatial trajectory is visualized. 

The perspective main view shows the projection of the actual trajectory on the 
orbital plane, rather than the trajectory itself. The reason for this is twofold. 
The orbital trajectory (with its typical cycloidal shape), when shown without lines 
projected on the orbital reference plane, is ambiguous and might seem to bend out of 
the orbital plane. This illusion results from the viewer's familiarity with common 
objects such as a coil spring, and has first been reported in reference 7. There- 
fore, the trajectory cannot be shown without its projection on the orbital plane. 
Second, the symbolic enhancements and burn vectors relate to the in-plane motion and 
match with the trajectory projection on the orbital plane. Thus, both the trajec- 
tory and its projection should be visualized. However, in a perspective plan view, 
i.e., viewed along the y° axis, both the trajectory and its projection on the 
orbital plane will show up as separate curves, which might be highly confusing. 
Therefore a compromise has been sought, in which the projection is shown together 
with pedestals placed at the way points and orthogonal to the orbital plane. The 
end of the pedestal shows the actual three-dimensional position of the way point. 

In spite of these restrictions, the proposed display clearly demonstrates the 
usefulness of interactive graphical trajectory design. The use of the graphical, 
symbolical, and computational enhancements indicates the direction in which a solu- 
tion for a multivehicle environment display should be sought. A still unanswered 
question relates to the degree of automation that should be introduced into the 
display. Parts of the mission could be performed through the use of optimization 
techniques; e.g., to find the fuel-optimal way point which clears a spatial 
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constraint in part of the mission, or to find a way point which satisfies the 
terminal constraints. However, since the solution space of a complex situation is 
virtually infinite, it is doubtful whether this mission can be performed entirely 
automatically. It is therefore expected that frequently occurring routine opera- 
tions, such as searching the local solution space for the optimal location of a way 
point, might be handed over to an optimization scheme. These solutions can be 
reviewed by the operator and manually changed if necessary. 

An experimental program is planned in which operators will carry out a series 
of design missions which vary in complexity and constraints. In a tutorial session, 
the operators will be familiarized with the orbital motions, orbital control meth- 
ods, and operational constraints. Furthermore, the operators will be familiarized 
with the system control functions of the viewing system motions and way point edit- 
ing process. Each operator action, i.e., viewing system or trajectory design, will 
be time-marked and recorded. Statistics of the viewing system actions will show 
preferred viewing situations for each condition. Review of the trajectory design 
actions might reveal the existence of heuristic design rules which might be utilized 
in automated design schemes. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, CA 94035, April 14, 1988 


APPENDIX 


MATHEMATICAL BACKGROUND 


Relative Motion in Space Station Coordinates 

The relative position and velocity of a co-orbiting spacecraft in Space Station 
coordinates are x° = {x 0 ,y°,z°} and x° = {x°,y 0 ,z°}, and are obtained from Space 
Station-based measurement equipment. Since the displacements and velocities in the 
out-of-orbital-plane direction y° are usually much smaller than the ones in the 
orbital plane in the x° and z° directions, and since an out-of-plane maneuvering 
burn does not significantly alter the total orbital energy, the out-of-plane motion 
can be decoupled from the in-plane motion. Hence, the in-plane motion will be 
analyzed first. 

The in-plane orbital situation is shown in figure 16. At the initial time 
t = t Q the radius of the Space Station orbit R Q is given by 

R = R_ + h ( A 1 ) 

o E 

where R E = 6,378,140 m is the equatorial Earth radius and h = 480,000 m is the 
altitude of the Space Station orbit above the Earth surface. The absolute orbital 
velocity of the Space Station is then given by 


V 

o 



(A2) 


where GM = 3.986005*10^ m^/sec 2 is the geocentric gravitational constant, 
simplicity the curvature of the V-bar is assumed to be negligible, so that 


For 


s » x 


r s -z 


(A3a) 

(A3b) 
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a 
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r , * ,2 ,*. 2 , 1/2 

[(s) + (r) ] 

tan~\r/s) 

tan~^(r/s) + 180° 


s > 0 
s < 0 


(A4a) 

(A4b) 


where s is the distance measured along the V-bar between the Space Station and the 
spacecraft's R-bar; r is the distance of the spacecraft above the V-bar, measured 
along the R-bar; v is the magnitude of the relative velocity; and a is the angle 
of its direction, measured from the V-bar in an upward direction (positive rotation 
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in the right-hand system). The relative velocity reflects the situation just after 
the activation of a maneuvering burn. 

The squared absolute orbital velocity follows from the geometry of figure 16 
and is then given by 


2 2 
V = r 
o 




The orbital energy equation can be written as 


2 

1 2 r 

a = R " GM 


(A5) 


(A6) 


where a is the semimajor axis of the spacecraft orbit. Substituting (A5) into 
(A6) and using (A2) yields 


(!)•»- [(t) 3 * a (U’fc) * * (fc) 


and using the fact that R = R 0 + r yields 


= 1 - 2 
a 


Note that in equation (A7b) the ratio r/R Q rather than R/R Q is used, 
done for reasons of numerical accuracy. 


(A7a) 


(A7b) 


This is 


The shape of the orbit is given by 


R = 


a( 1 - e ; 

1 + e cos v 


(A8) 


where v is the true anomaly, which is the angular position in orbit after perigee 
passage and e the eccentricity of the orbit. Defining the auxiliary angle E, the 
eccentric anomaly, and manipulating the orbital energy equation and the orbital 
angular momentum equation yields (see ref. 1) 
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(A9a) 

(A9b) 


Now, for known v, a, r, V Q , and R Q , the values for A and B can be computed by 
first evaluating R/a with equation (A7b) and substituting this value in equations 
(A9a,b). The eccentricity follows from 


e 



♦ B 2 ] 1/2 


(A10) 


and E from 


E = tan -1 (B/A) 


- 1 , 


A > 0 

E = tan” 1 (B/A) + 180° A < 0 
The mean anomaly at time t = t Q is computed from Kepler's equation 

M(t Q ) = n(t Q - T) = E - e sin E = E - B 


(All) 


(A12) 


where n is the mean motion or angular orbital rate of the spacecraft in radians 
per second, given by 


n 



(A13) 


and t = T is the time of perigee passage. The mean motion n Q of the Space 
Station is given by 


n 


o 



Combining equations (A13) and ( A14) yields 


n = n C 
o 


3/2 


where 


r (R/a) 

" (R/R Q ) 


(A14) 


(A15a) 


(A15b) 


and the difference in mean motion An = n - n Q between spacecraft and Space Station 
is approximated by 
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Figure 5.- Orbital motion after V-bar burn, (a) Shape of the orbit 
(b) Trajectory relative to the Space Station. 
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Figure 6.- Example of a three-burn maneuver. 
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Figure 4.- Orbital motion after R-bar burn, (a) Shape of the orbit 
(b) Trajectory relative to the Space Station. 
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- Orbital motion for an initially stationary spacecraft with radial separa- 
(a) Shape of the orbit, (b) Trajectory relative to the Space Station. 
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Figure 1.- Definition of the orbit, (a) Orientation of the orbital plane, 
(b) Shape of the orbit in the orbital plane. 
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1/2 


F = 


2 "< fc 0 >‘ 
“ (t °, * ~f 


<fr = tan -1 tnw(t 0 )/w(t 0 )] 

4» = tan _1 [nw(t )/w(t )] + 180° 
o o 


w(t Q ) > o 
w(t Q ) < 0 


(A48a) 

(A48b) 


Using equation (A47), equations (A45a,b) can be reformulated to include the angle <t> 

w(t) = F sin[n(t - t Q ) + ♦] (A49a) 

w(t) = Fn cos[n(t - t Q ) + 4>] (A49b) 

The relative out-of-plane motion at t = t is then computed by substituting the 
values of F and <t> in equations (A49a,b). 


Inverse Computations for Out-of-Plane Motion 

The inverse procedure for out-of-plane motion has a simple closed-form solu- 
tion. As is the case for in-plane motion, the problem is to find the value of 
w(t Q ) which, starting from initial position w(t Q ), results in a final position 
w(t) at t = t. Equation (A49a) is written in expanded form 

w(t) = F cos 4> sin n(t - t ) + F sin <t» cos n(t - t ) (A50) 

o o 

Substituting equations (A46a,b) into equation (A50) and rearranging them yields the 
required w(t Q ) 

w(t ) = n[w(t) - w(t )cos n(t - t )]/sin n(t - t ) (A51) 

o o o o 
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R q which are defined by the initial relative position at t = t Q . The correspond- 
ing orbital parameters and relative position are computed with equations (A7-A21). 

Next, using equations (A44a,b) the derivatives of R/a are computed by using 
the same values of v and o of the previous step as were used in the forward algo- 
rithm. The derivatives of R/a are used in equations (42a, b) and (43a, b) to com- 
pute the derivatives of A and B. The derivatives of a, n, An, e, E Q , and M Q are 
computed with equations (A33-A41). Then, together with the values of D 1 and D 2 
computed with equations (A28a,b), all computed derivatives are used to compute the 
derivatives of and D 2 with equations (A32a,b). Finally, all computed values 
are used in equations (A29a,b) to compute the derivatives of s and v. 

The 2x2 Jacobian matrix can now be formed. The gradient vector of equa- 
tion (A24) is computed by simple matrix operations. In the computation of the 
second-derivative matrix, the second-order term on the right-hand side of equa- 
tion (A25) can be neglected. The new value of u is computed with equation (A23). 


Out-of-Plane Motion 

The out-of-plane motion of a co-orbiting spacecraft results from a noncoplanar- 
ity between the orbital plane of the spacecraft and the plane of the Space Sta- 
tion. An out-of-plane maneuvering burn alters the relative inclination of the 
orbital plane of the spacecraft and thus the out-of-plane motion. The out-of-plane 
displacement w and velocity w at time t are described by 

w(t) = y° = F sin n(t - t) (A45a) 

w(t) = y° = Fn cos n(t - t) (A45b) 

where F is the maximum out-of-plane deviation from the Space Station orbital 
plane, t = t is the time at which the Space Station orbital plane is crossed in the 
positive y° direction, and n is the mean motion of the orbit. 

At the initial time t = t Q the initial out-of-plane position and velocity are 
given by 


w(t ) = F sin $ (A46a) 

o 

w(t Q ) = Fn cos <f> (A46b) 

where $ is the orbital angle after plane crossing at t = t Q defined as 

<D = n(t Q - t) (A47) 

and w(t Q ) reflects the velocity just after a maneuvering burn at t = t Q . For a 
known initial position and velocity, F and 4> can be computed using equa- 
tions (A46a,b) 
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Differentiating equation (A40) yields 


3a 

3x 



(A41) 


Derivatives of A and B 


The derivatives of A and B will be evaluated for v and a separately. 
Differentiating equation (A9a) yields 

3A _ 3(R/a) 

3V 3v 


3A 3(R/a) 

3a ~ 3a 



(A42a) 

(A42b) 

(A43a) 

(A43b) 


Derivatives of R/a 

The derivatives of R/a are obtained by differentiating equation (A7). 

(A44a) 
(A44b) 


3(R/a) 

3V 

3 (R/a) 
3a 


■ l (0 ™ ‘ * ( 0(0 

4-J(r) • 


Implementation of the Inverse Algorithm 

The inverse algorithm is easily implemented by computing the derivatives and 
dependencies in the reversed order as presented here. Computer coding can be highly 
simplified by identifying compound terms which appear more than once in the pro- 
cess. These terms have to be computed only once and their values are substituted 
whenever necessary. The above equations are presented in a format which makes 
identification of these terms easier. 

The algorithm starts with the computations of the forward algorithm, i.e., 
based on the values of v and a of the previous step and on the values of R and 
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The t derivative of E Q involves the derivative of an arc tangent. The general 
form for this derivative is given by 


3(tan~^z) 

ax 


1 3z 

. 2 3x 

1 + z 


(A36) 


Substituting B/A for z in equation (A36) and using equation (A10) yields 


3E 

o 


3X 




(A37) 


Thus, evaluating the derivative of M Q also requires the derivatives of A and B. 
These derivatives will be evaluated later. 


Derivatives of e 


The expression for the derivatives of e is obtained by differentiating equa- 
tion (A10) 


|s = i A fi& + b isl 

3x e lax 3xJ 


(A38) 


Derivatives of n and An 


derivatives of 
the term R/R, 


The expression for the 
equation (A15a). Note that Vll\y WVi 111 11/ IIq 

R = R q + r is the radius of the orbit at 
term R/a 
a. Since 
given by 


n and An is obtained 
does not depend on v 


by differentiating 
or a, since 


t = t 


o* 


This is in contrast with 


in which a (which determines the shape of the orbit) depends on 


the 
v and 


n. 


is invariant, the derivatives of n and An are equivalent and are 


3n 

3x 


3 An 

3X 


= 1.5n 


[^]/@ 


(A39) 


The derivatives of the term R/a will be evaluated last. 


Derivatives of a 


a is written as 



(AH0) 
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and accordingly 


3 sin M.j 
ax 




cos M 1 


(A31a) 


a sin M 
o 

ax 


3M 

t “ COS M 
ax o 


(A3 lb) 


a cos M.j 
ax 




sin M 1 


a cos M 
o 

ax 



sin M 

o 


(A31c) 
(A3 Id) 


Differentiating equations (A28a,b), substituting equations (A31a-d), and using 
equations (A28a,b) again yields 


a D 1 
ax 


3M 

ax 


2 °2 + If (t - fc o )cos M 1 


(A32a) 


a d 5 aM 

it - sr D i ♦ If (t - ‘o' 51 " M i 


(A32b) 


Now, in order to evaluate the derivatives of equations (A29-A32), the derivatives 
of M Q , e, n, An, and a are required. 


Derivatives of M Q 

M q follows from equation (A12) and is given by 

M = E - B (A33) 

o o 

where E Q is the eccentric anomaly at time t = t Q . Differentiating equation (A33) 
yields 


Since 



3x ax 3x 


(A34) 


E = tan -1 (B/A) 
o 


(A35) 
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(A26) 


dS 

3s" 

3V 

3a 

3 r 

3r 

,3v 

3a_ 


The sequence of dependencies for computing the derivatives of the Jacobian matrix of 
equation (A26) are outlined here. 


Derivatives of s and r 


Equations (A19a,b) of the relative trajectory can be written as 



s(t) = s(t Q ) + R Q [An(t - t Q ) + 2eD^] 

(A27a) 

where 

r(t) = r(t ) + aeD_ 
0 2 

(A27b) 


D 1 = sin M 1 - sin M q 

( A28a ) 


D,, = -cos M, + cos Nf 
2 1 o 

(A28b) 


It is clear that each one of the terms on the right-hand side of equations (A27a,b) 
depends on both v and a. Accordingly, the derivatives of these elements with 
respect to v and a have to be computed. Differentiating equations (A27a,b) yields 


If = "oj^ 1 - V + 2 If D 1 + 2e ^ir] U29a) 

If ■ [ft * * a lf] D 2 * ae (A29b) 

where the symbol x is used to denote either v or a. 

Derivatives of D 1 and D 2 

Following equations (A21a,b), can be written as 

M. = n(t - T) = n(t - T) + n(t - t ) = M + n(t - t ) (A30) 

I o o o o 
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where 


M 1 = n(t - T) 

«o 2 n(t o - T) 


and T is obtained from equation (A12). 


(A21a) 

(A21b) 


The Inverse Algorithm 


It has been shown that in order to give the operator direct control over the 
position of way points, the relative trajectory equations have to be solved back- 
ward; i.e., for a given initial relative location at t = t Q and a final location 
at t = t, the corresponding values of v and a at t = t Q have to be found. A 
Newton-Raphson scheme has been chosen to solve this problem. 

Define the state of the system at t = t by the vector x = {s,r} and the 
control applied at t = t Q to bring the system to this state by the vector 
u = {v,a}. Then equations (A1)-(A21) show that x is a sophisticated nonlinear 
function of u. The problem is to find u, which brings x to a given reference 
position x r . ^ dismatch score is defined as 

J(u) = (x - x r ) T (x - x r ) (A22) 


For the 


ith iteration 


Si 


is computed from u^_i according to 



3 2 J 

-1 

3 J 

1 

3U 2 

-i-1 

8u 


(A23) 


The gradient vector is given by 


3J 

3U 


*1-1 


T 9 21 

= (x - x )* — 
- -r 3u 


(A24) 


^i-1 


and the 2x2 second-derivative matrix is given by 

vT 

1 ,n 

2 


3 2 J 


3U 


( 3X 3x 

iu/ iu 


iii-1 


T a2 - 

+ (21 - 2lJ T ~ 2 

Sii-1 


(A25) 


Si-1 


Both the gradient vector and the second-derivative matrix require the computation of 
the 2x2 Jacobian matrix, which is given by 
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An » -1 .5n( 1 - C) 


(A16) 


| 

i 


i 


i 


i 


Note that for numerical reasons it is preferred to use a first-order Taylor series 
approximation in the computation of An. 

It is clear from the above derivations that for a given initial relative posi- 
tion and velocity at t = t Q , the orbital parameter set, which characterizes the 
shape of the orbit, can be fully determined. This parameter set can now be used to 
compute the relative position and velocity at time t = t. This is done in the 
following steps. First, the mean anomaly at time t is computed from 

M(t) = M(t Q ) + n(t - t ) ( A 1 7 ) 

Now, Kepler's equation (eq. (12)) can again be employed to compute the eccentric 
anomaly at time t = t, E(t). However, a simple closed-form solution for solving 
E does not exist. The common practice in classical celestial mechanics is to use 
Fourier-Bessel functions and series expansions (see ref. 1). An alternative is to 
solve the equation numerically by an iterative process such as Newton's false-root 
method. Once the eccentric anomally E(t) is computed, the ratio R/a at t = t 
can be found with equation (A9a). Multiplying the ratio with the known value for 
a and subtracting R Q yields the relative distance r. By using equation (A8) the 
true anomaly v(t) at t = t can be found. The relative distance s is obtained 
from 

s(t) = R Q [v(t) - v Q (t)] (A18) 

where v Q is the true anomaly of the Space Station. 

In this research, the Fourier-Bessel series expansions are used (see ref. 1). 
Since the eccentricity of the orbit is very small, second- and higher-order terms 
were dropped. The resulting equations for the relative trajectories and trajectory 
rates, starting from t = t Q are given by 

s(t) = s(t Q ) + R Q [ An(t - t Q ) + 2e(sin M 1 - sin M q )] (A19a) 

r(t) = r(t ) + ae(cos M - cos M.) (A19b) 

O 0 1 

and 

s(t) = R Q [An + 2en cos M^] (A20a) 

r(t) - aen sin M 1 (A20b) 


> 
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Figure 11 



Time-axis format with operational constraints and day-night indication 



Figure 12.- Departure constraints. 



Figure 13.- Arrival constraints. 



Figure 14.- Plume impingement constraints. 
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